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Abstract. The Harmonic Balance method provides a heuristic approach for 
finding truncated Fourier series as an approximation to the periodic solutions 
of ordinary differential equations. Another natural way for obtaining these type 
of approximations consists in applying numerical methods. In this paper we 
■ recover the pioneering results of Stokes and Urabe that provide a theoretical 

O^l , basis for proving that near these truncated series, whatever is the way they 

• have been obtained, there are actual periodic solutions of the equation. We will 

O , restrict our attention to one-dimensional non-autonomous ordinary differential 

equations and we apply the results obtained to a couple of concrete examples 
00 ■ coming from planar autonomous systems. 
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1. Introduction and main results 

Consider the real non-autonomous differential equation 

x' = X{x,t), (1) 

where the prime denotes the derivative with respect to t, X : Q x [0, 27r] — t- M is a 
C^-function, 27r-periodic in t, and C M is a given open interval. 

There are several methods for finding approximations to the periodic solutions 
of ([I]). For instance, the Harmonic Balance method (HBM), recalled in subsec- 
tion [2TT1 or simply the numerical approximations of the solutions of the differential 
psj ■ equations. In any case, from all the methods we can get a truncated Fourier se- 

O . ries, namely a trigonometric polynomial, that "approximates" an actual periodic 

solution of the equation. The aim of this work is to recover some old results of 
Stokes and Urabe that allow to use these approximations to prove that near them 
there are actual periodic solutions and also provide explicit bounds, in the infinity 
norm, of the distance between both functions. To the best of our knowledge these 
i results are rarely used in the papers dealing with HBM. 

When the methods are applied to concrete examples one has to deal with the 
coefficients of the truncated Fourier series that are rational numbers (once some 
number of significative digits is fixed, see the examples of Section |4]) that make 
more difficult the subsequent computations. At this point we introduce in this 
setting a classical tool, that as far as we know has never been used in this type 
of problems: we approximate all the coefficients of the truncated Fourier series by 
suitable convergents of their respective expansions in continuous fractions. This is 
done in such a way that using these new coefficients we obtain a new approximate 
solution that is essentially at the same distance to the actual solution that the 
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starting approximation. With this method we obtain trigonometric polynomials 
with nice rational coefficients that approximate the periodic solutions. 

Before stating our main result, and following 0,15], we introduce some concepts. 
Let x{t) be a 27r-periodic C^-function, we will say that x{t) is noncritical with 
respect to ([1]) if 

j^^ ^X{xit),t)dty^O. (2) 

Notice that if x{t) is a periodic solution of ([T]) then the concept of noncritical is 
equivalent to the one of being hyperbolic, see [3] . 

As we will see in Lemma [2.11 if x{t) is noncritical w.r.t. equation ([1]), the linear 
periodic system 

y' = ^X{x{t),t)y + b{t), 

has a unique periodic solution yb{t) for each smooth 27r-periodic function b{t). 
Moreover, once X and x are fixed, there exists a constant M such that 

H\L < M\\b\\2, (3) 

where as usual, for a continuous 27r-periodic function /, 




'^Pm, ||/|L=max|/(x)| and ||/||, < 



Any constant satisfying ([3]) will be called a deformation constant associated to x 
and X. Finally, consider 

s{t) := x{t) - X{x{t),t). (4) 

We will say that x(t) is an approximate solution of ([T]) with accuracy S = \\s\\2. 
For simplicity, if S* > S* we also will say that x{t) has accuracy S. Notice that 
actual periodic solutions of ([T]) have accuracy 0, in this sense, the function s{t) 
measures how far is x{t) of being an actual periodic solution of ([1]). 

Next theorem improves some of the results of Stokes [3] and Urabe [B] in the 
one-dimensional setting. More concretely, in those papers they prove the existence 
and uniqueness of the periodic orbit when AM'^KS < 1. We present a similar proof 
with the small improvement 2M'^KS < 1. Moreover our result gives, under an 
additional condition, the hyperbolicity of the periodic orbit. 

Theorem 1.1. Let x(t) be a 271 -periodic -function such that: 

- it is noncritical w.r.t. equation ([T]) and has M as a deformation constant, 

- it has accuracy S w.r.t. equation ([T]). 

Given I := [min^t^^-j x(t) — 2MS', maxj^gK} + 2MS] C Q, let K < oo be a 
constant such that 

^\x{x,t) 



max 

(x,t)e/x[0,27r] 



dx'^ 



< K. 



Then, if 

2M'^KS < 1, 

there exists a 211 -periodic solution x*{t) of ([T]) satisfying 

\\x*-x\\ <2MS 
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and it is the unique periodic solution of the equation entirely contained in this 
strip. If in addition, 

'2% 



r2n Q 



then the periodic orbit x*(t) is hyperbolic and its stability is given by the sign of 
this integral. 



Once some approximate solution is guessed, for applying Theorem ll.il we need 
to compute the three constants appearing in its statement. In general, K and 5* 
can be easily obtained. Recall for instance that ||s||2, when s is a trigonometric 
polynomial, can be computed from Parseval's Theorem. On the other hand M 
is much more difficult to be estimated. In Lemma 12.31 we give a result useful for 
computing it in concrete cases, that is different from the approach used in [SllHlE]. 

Assuming that a non-autonomous differential equation has an hyperbolic peri- 
odic orbit, the results of [B] also guarantee that, if take a suitable trigonometric 
polynomial f{t) of sufficiently high degree, we can apply the first part of The- 
orem 11.11 Intuitively, while the value of the accuracy S goes to zero when we 
increase the degree of the trigonometric polynomial, the values M and K remain 
bounded. Thus at some moment it holds that 2M'^KS < 1. 

In Section |4] we apply Theorem 11.11 to localize the limit cycles and prove its 
uniqueness, in a given region, and its hyperbolicity for two planar polynomial 
autonomous systems. The first one is considered in Subsection 14.11 and is a simple 
example for which the exact limit cycle is already known. We do our study step by 
step to illustrate how the method suggested by Theorem 11.11 works in a concrete 
example. In particular we obtain an approximation x{t) of the periodic orbit by 
using a combination between the HBM until order 10 and a suitable choice of 
the convergents obtained from the theory of continuous fractions applied to the 
approach obtained by the HBM. 

The second case corresponds to the rigid cubic system 

X = -y + ^(1 - X - lOx^), 
y = X + ^{1 - X - lOx^), 

that in polar coordinates writes as f = r/10 — cos(^)r^/10 — cos^(^)r'^, ^ = 1, or 
equivalently, 

which has a unique positive periodic orbit, see also [2]. Notice that we have 
renamed 6 as t. We prove: 

Proposition 1.2. Consider the periodic function 

m = ^ - ^ cos(t) - 1 sin(t) - ^ cos(2i) - 1 sin(2t) - ^ cos(3t). 
Then, the differential equation (|5]) has a periodic solution r*{t), such that 

||r < 0.042, 

which is hyperbolic and stable and it is the only periodic solution of (|5]) contained 
in this strip. 
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As we will see, in this case we will find computational difficulties to obtain the 
order three approximation given by the HBM. So we will get it first approaching 
numerically the periodic solution; then computing, also numerically, the first terms 
of its Fourier series and finally using again the continuous fractions approach to 
simplify the values appearing in our computations. We also will see that the same 
approach works for other concrete rigid systems. 

Similar examples for second order differential equations have also been studied 
in [3. 



2. Preliminary results 



This section contains some technical lemmas that are useful for proving The- 
orem 11.11 and for obtaining in concrete examples the constants appearing in its 
statement. We also include a very short overview of the HBM adapted to our 
interests. See for a more general point of view on the HBM. 

As usual, given A C M, 1a : — denotes the characteristic function of A, 
that is, the function takes the value 1 when x & A and the value otherwise. 

Lemma 2.1. Let a{t) andhit) he continuous real 2ti -periodic functions. Consider 
the non- autonomous linear ordinary differential equation 



x' = a{t)x + hit). 



(6) 



If A(2n) 7^ 0, where A{t) := Jl^a{s)ds, then for each h{t) the equation ([6]) has a 

unique 2tx -periodic solution Xb{t) := J^^ H(t, s)h{s)ds, where the kernel H(t,s) is 
given hy the piecewise function 



H{t,s) 



oA{t) 



1 _ eA(27r) 

Moreover \ \xb\\^ < 27rmaxtg [0,277] \ \H{t 
Proof. Since is linear, its general solution is 



[e-^(^)l[o,](s) + e^(^-)-^(^)l[,2.](3)] 



(7) 



x{t) 



Xo+ h{s)e-^^'^di 



If we impose that the solution is 27r-periodic, i.e., x(0) = x{2tt), we get 



,A{27r) 



Xo 



2lT 



h{s)e-^^'Us. 







1 _ e^(2T) 

By replacing xq in ([8]) by the right hand side of ([9]) we obtain that 

eA(t) r r2n 



(8) 



(9) 



Xb{t) 



1 _ (,A(2n) 



1 _ eA(27r) 
2tt 

H{t,s)h{s)ds. 




Jo Jo 

gA(2.) r h(^s)e-^^'^ds+ f h{s)e-^^'^ds 
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Therefore the first assertion follows. On another hand, by the Cauchy-Schwarz 
inequality, 

/ j-2lT / j-2-K 

\xb{t)\ <J J H^it,s)dsJ J b^{s)ds. 

Therefore 

IkbIL <27r max ||if(t,-)|U|6||,. 
This complete the proof. □ 
Corollary 2.2. A deformation constant M associated to x and X is 

M:=27i max 

t6[0,27r] 

where H is given in ([7]) with A{t) = ■^X(x{t),t) dt. 

Now we prove a technical result that will allow us to compute in practice defor- 
mation constants. In fact we will find an upper bound of M that will avoid the 
integration step needed in the computation of the norm || ■ ||2- First, we introduce 
some notation. 

Given a function A : [0, 27r] — > M, a partition tj = ih,i = 0,1, . . . , N, of the 
interval [0, 27r], where h = 27t/N, and a positive number i, we consider the function 
L : [0, 27i] — )■ M given by the continuous linear piecewise function joining the points 
{ti,A{ti) -i). Notice that L{t) = ^iJo^ Li{t)li^, where It = and 

We will say that L is an adequate lower bound of A if it holds that L{t) < A{t) 
for all t G [0, 27r]. It is clear that smooth functions have always adequate functions, 
that approach to them. 

In next result we will use the following functions 

■m-l 7V-1 a 

gPrr- 



^Ut) := 5^ + Yl ^' + (1 - (e"""* - e"'"*'") , (10) 

1=0 i=m—l 

where 

Ju Ju «i 



and A = e^(2-). 

Lemma 2.3. Let L he an adequate lower hound of A, where A is the function 
given in Lemma \2.1\ Consider the functions '^m{t), m = 0,1, . . . , N — 1, given 
in (fTOj) . Therefore, following also the notation introduced in that Lemma, it holds 
that I 1^ < I I2, where 



N = ^^ max e^W, 



— — iiia-vv o A 

|1 — A| i6[0,27r] \ 



m=0 
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Proof. Recall that from Lemma \2.1\ \\xi,\\^ < M \ \b\\^, where 

M:=2n m8.x \\H{t, ■)\\,. 

t6[0,27r] 



So we will find an upper bound of M. Since 

1 _ eA(27r) 



it holds that 

1 , 

yZn I i — A| 

where 

G{t) := Te-^^Wrfs + A^ r e-^^^^^ds< f e-'^^^'Us + T^-^^^^^^is, 
Jo Jt Jo Jt 

because L{t) < A{t), for all t e [0,27r]. 

Assume that t & I^- Then 



•^^ i=m i=m-l "^*™ 

Therefore, for t E Im, 



e 



G(t)<5^J, + A2 ^ j, + (1-A2) / e^-^+'^-rfs = vl/„(t). 

j=0 i=m~l 

As a consequence, for t G [0, 27r], 

G(t)<5^vi/„(t)i,^(t), 

m=0 

and the result follows. □ 

Remark 2.4. Notice that the above lemma provides a way for computing a defor- 
mation constant where there is no need of computing integrals. This will be very 
useful in concrete application, where the primitive of e'"^^^*^ is not computable and 
so Corollary \2.i^ is difficult to apply for obtaining M. 



In next result, which introduces the constant K appearing in Theorem II ■![ D° 
denotes the topological interior of D. 

Lemma 2.5. Consider X as in ([T]). Let D be a closed interval and let x{t) be a 
27r -periodic -function, such that {x{t) : t G M} C D° . Define 

d 

R{z, t) := X{x{t) + z,t)- X{x{t),t) - —X{x{t),t)z, (11) 

ox 

for all z such that {x{t) + 2 : t G M} C D. Then 
{{} \Riz,t)\<f\z\\ 
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(a) \R{z,t) — R{z,t)\ < K max{\z\, \z\) \z — z\, 
where 



K 



max 

(x,t)eDx[0,27r] 



X{x,t) 



Proof, (i). By using the Taylor's formula, for each t it holds that 

Xixit) + z,t)= X{x{t),t) + ^X{x{t),t)z + ]^^X{i{t),t)z' 

for some E,{t) G {x{t),x(t) + z). Therefore 

1 92 



\R{z,t)\ 



2 dx' 



\z\'< 



K . 



as we wanted to prove. 

{ii). From Rolle's Theorem for each fixed t it follows that there exists ri{t) e (2;, z) 
such that 



\R{z,t) - Riz,t)\ < 



d 



dz 

dz 



Rivit),t) 



\z — z\. 



Applying again this theorem, but now to -§^R, noticing that ^R{z^ t) 
obtain that 

d 



z=0 



0, we 



dz 



R{vit),t) 



< 



dz'- 



■R{uj{t),t) 



dx' 



:X{u;{t),t) 



m\<Km\, 



where u(t) e (0,?7(t)). Note also that 

|?7(t)| < max(|2;|, l^]). 
Hence, the result follows combining the three inequalities. 



□ 



2.1. The Harmonic Balance method. In this subsection we recall the HBM 
adapted to the setting of one-dimensional 27r-periodic non-autonomous differential 
equations. 

We are interested in finding periodic solutions of the 27r-periodic differential 
equation ([T]), or equivalently, periodic functions which satisfy the following func- 
tional equation 

J^{x{t)) := x'{t) - X{x{t),t) = 0. (12) 

Recall that any smooth 27r-periodic function x{t) can be written as its Fourier 
series, 



x{t) = + (amcos(mt) + bjnSm{mt)) 



m=l 



where 



2lT 



x(t) cos(mt) dt, and b„ 



TT 



2tt 



x{t) sin(mt) dt, 



for all m > 0. Hence it is natural to try to approach the periodic solutions of 
the functional equation f|T2l) by using truncated Fourier series, i.e. trigonometric 
polynomials. 
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Let us describe the HBM of order A^. Consider a trigonometric polynomial 

N 



y_ 



(t) = ^ + ^ {vm cos(mt) + Sm sin(mt)) 



771=1 



with unknowns = rm{N),Sm = Sm^N) for all m < N. Then compute the 
27r-periodic function J^{yj^{t)). It has also an associated Fourier series 

J^iVNi^)) = ^ + i-^rn cos(mt) + Br,^ sin(mt)) , 

m=l 

where Am = Am{'r,s) and Bm = Bm{'r,s), m > 0, with r = (ro, ri, . . . , r^) and 
s = (si, . . . , s^). The HBM consists in finding values r and s such that 

^^(r, s) = and B,n{r, s) = for 0<m<N. (13) 

The above set of equations is usually a very difficult non-linear system of equa- 
tions and for this reason in many works, see for instance |1] and the references 
therein, only small values of are considered. We also remark that in general the 
coefficients of and do not coincide at all. 

Notice that equations f lT5]l are equivalent to 

/'27r /'2n 

/ T{y^{t)) cos(mt) dt = Q and / J^{yj^{t)) sin(mt) dt = 0, 
Jo Jo 

for < m < iV. 

The hope of the method is that the trigonometric polynomials found using this 
approach are "near" actual periodic solutions of the differential equation ([T]). In 
any case, as far as we know, the BHM for small is only a heuristic method that 
sometimes works quite well. 

To end this subsection, we want to comment a main difference between the non- 
autonomous case treated here and the autonomous one. In this second situation 
the periods of the searched periodic orbits, or equivalently their frequencies, are 
also treated as unknowns. Then the methods works similarly, see again [1]. 

3. Proof of the main result 



Proof of Theorem \l.l[ As a first step we prove the following result: consider the 
nonlinear differential equation 

z = X{z + x{t),t) - X{x{t),t) - s{t), (14) 

where s{t) is given in (j4]). Then a 27r-periodic function z{t) is a solution of (fT4|) if 
and only if z{t) + x{t) is a 27r-periodic solution of ([1]). 
This is a consequence of the following equalities 

{z{t) + x{t)y =[X{z{t) + x{t),t) - X{x{t),t) - s{t)] + [X(x(t), t) + s{t)] 
=X{z{t) + x{t),t). 

By using the function 

d 

R(z, t) = X(z + x(t),t) - X(x(t),t) - —X(x(t),t)z, 

ox 
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introduced in Lemma I2.5[ equation f|T4l) can be written as 

d 

z' = —X(x(t),t)z + R(z,t)-s(t). (15) 
ox 

Let V be the space of 27r-periodic C^-functions. To prove the first part of the 
theorem it suffices to see that equation f lT^ has a unique C^, 27r-periodic solution 
z*{t), which belongs to the set 

^={zeV: II^IL < ^MS}. 

To prove this last assertion we will construct a contractive map T : A/" — )■ A/". 
Because A/" is a complete space with the || ■ ||^ norm, its fixed point will be a 
continuous function in A/" that will satisfy an integral equation, equivalent to (|T5l) . 
Finally we will see that this fixed point is in fact a function and that it satisfies 
equation (|T5i) . 

Let us define T. If z E J\f then T{z) is defined as the unique 27r-periodic solution 
of the linear differential equation 

y' = ^X{x{t),t)y + R{z{t),t) - s{t). 

Notice that this map is well defined, by Lemma 12.1^ because x{t) is noncritical 
w.r.t. equation ([T]). Then zi satisfies 

d 

z[ = —X{x{t),t)z^ + R{z{t),t) - s{t). 

Let us prove that T maps A/" into A/" and that it is a contraction. By Lemmas 
12.11 and 12.51 and the hypotheses of the theorem 

\\nz)\L = IkilL < M\\R{z{-), ■) - s{-)\U <M{\\R{z{.), + S) 

<M{\\R{z{.),.)\L + S)<M{^\\z\\l + S) 

< M{2KM'^S^ + S) < 2MS, 

where we have used in the last inequality that 2M^KS < 1. 

To show that T is a contraction on JV, take z,z E M and denote by zi = T{z), 
z\ = T{z). Then 

z[ = ^X{xit),t)z, + R{z{t),t) - sit), 
d 

z[ = —X{x{t),t)zi + R{z{t),t) - s{t). 



Therefore 



{zi - = ^X(x(t), - + Riz{t),t) - R{z{t),t). 



Again by Lemmas 12.11 and 12.51 and the hypotheses of the theorem, 

\\T{z) - T{-z)\L =lki - ZilL < M\\Riz{-), ■) - R{z{.), OIL 

<Mi^'max(||^|L, II^IDIIz - z\L < 2M'KS\\z - z\ 

as we wanted to prove, because recall that 2M'^KS < 1. 
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Therefore the sequence of functions {zn(t)} defined as 

4+1 W = ^X{x{t),t)Zn+l{t) + R{Zn{t),t) - S{t), 

with any Zgit) G TV, and Zn+i(t) chosen to be periodic, converges uniformly to 
some function x* [t) G J\f. In fact we also have that 



Zn+lit) = Zn+l{0) + 

Therefore 

x*{t) = x*{0) + 



d_ 

dx 



X[x{w),w)Zn+l{w) + R{Zn{w),w) — s{w) ) dw . 



X {x{w) , w)x* (w) + R{x*{w),w) — s{w) J dw. 



d_ 

dx 

We know that x*(t) is a continuous function, but from the above expression we 
obtain that it is indeed of class C^. Therefore x*{t) is a periodic solution of f llSp 
and is the only one in A^, as we wanted to see. 

To prove the hyperbolicity of x*(t) it suffices to show that 



—X{x*{t),t)dt^O, 



and study its sign, see [3]. We have that, fixed t, 



d 



d 



^^X{x*{t),t) = —X{x{t),t) + — - 

for some ^{t) G {x* {t) , x{t)) . Therefore, since we have already proved that \x*{t) — 
x{t)\ < 2MS, 



< 2KMS. 



Then 



27r 



dx 



X{x{t),t)dt 



2tt 



d_ 

dx 



X{x*{t),t)dt 



< inKMS < 



27r 

M 



and the results follows because by hypothesis the first integral is, in absolute value, 
bigger that 27i/M. □ 



4. Applications 

In this section we apply our result to prove the existence and localize a hyper- 
bolic limit cycle of some planar systems, which after some transformations can be 
converted into differential equations of the form ([T]). In the first case, although 
we know explicitly the limit cycle, we first use the HBM to approximate it and 
then Theorem 11.11 to prove in an alternative way its existence. In the second case 
we consider a planar rigid system. First, we found numerically an approximation 
of the limit cycle and from this approximation we propose a truncated Fourier 
series as a simpler approximation. Finally, Theorem II. II is used again to prove the 
existence and localize the limit cycle. 



2 , , .„.2^ (16) 
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4.1. A simple integrable case. Consider the planar ordinary differential equa- 
tion 

X = —y + x{a + dx"^ + exy + fy"^) 
y = X + y{a + dx"^ + exy + fy"^) 

In polar coordinates it writes as 

r = ar + {dcos^{6) + e sin(6') cos(6') + / sin2(6'))r^, ^ = 1, 
or equivalently, 

dv 

r' = — = ar + (rfcos^(t) + esin(t) cos(t) + /sin^(t))r^ := X(r,t), 

(JjL 

where we have renamed 6 as t. The above equation is a Bernoulli equation that 
can be solved explicitly. For simplicity we fix a = —1, d = 3, e = 2 and / = 1. 
Then we have the equation 

r = -r + (cos(2t) + sin(2t) + 2)r\ (17) 

Its solutions are r{t) = and 

1 

r{t) = : 



+ cos(2t) + A;e2* 

Therefore its unique positive periodic solution, which corresponds to the only limit 
cycle of f|T6l) for the given values of the parameters, is given by the ellipse 

r*(t) = , ^ (18) 
^/2T^^^ 

Moreover since 

/ —X{r*{t),t)dt = An>Q 
Jo or 

it is hyperbolic and unstable, see [3]. Its Fourier series is 

oo 

r*{t) = ^ + J2<^2kCos{2t), (19) 



where 



k=l 

4K 



ao = ^ ^ 1.491498374, ao/2 ^ 0.745749187, 
02 = i^^^ ^ -0.2016837219, 

V OTT 

^ -32g+20i^ _ 0.04065713288, 

V OTT 

as = 476g-296i^ _ -0.009092598292, 
_ -10624^+66047^ ^ 0.002133790322, 

V OTT 

^ i0M48^^6560M ~ _0.0005148662408, 

being = _R'(\/6/3) and E = E{\/6/3) the complete elliptic integrals of the first 
and second kind respectively, see pQ. 
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Let US forget that we know the exact solution and its full Fourier series to 
illustrate how to use the HBM and Theorem 11.11 for equation flTTl) to obtain an 
approach to the actual periodic solution flTSl) . 

Following the HBM, see subsection 12.1^ consider the equation 

jr(r(t)) = r'{t) + r{t) - (cos(2t) + sin(2t) + 2)r%t) = 0, (20) 

which is clearly equivalent to (fT7|) . 

Searching for a solution of the form r(t) = Tq and imposing that the first har- 
monic of J^{r{t)) vanishes we get that Tq + 2rQ = 0. The only positive solution of 
the equation is tq = \/2/2 ^ 0.7071 and this is the first order solution given by 
HBM. 

Motivated by the symmetries of f lT7|) for applying the second order HBM we 
search for an approximation of the form 

r{t) = Tq + r2 cos(2t). 

The vanishing of the coefficients of 1 and cos(2t) in the Fourier series of J^{r[t)) 
give the non-linear system: 

giro, ra) := tq - 2rl - ^rsrp - 3rlro - ^r^ = 0, 

9 3 

h{ro, ra) := r2 - - 6r2r^ - -r^ro - -r^ = 0. 

Doing the resultants Res((7, h, ro), Res(5', h, ra) we obtain that the solutions of the 
above system are also solutions of 

219720r^ - 18852r^ + 4269r^ - 328rl + 8 = 0, 

49437r^ - 70956r^ + 30708r^ - 4288r^ + 128 = 0. 

One of its solutions is ro ^ 0.7440456581 =: tq, ^ -0.2013905597 =: ^2- 

To know the accuracy of the periodic function f(t) = + cos(2t) as a solution 
of f lT7|) we compute 

S = I |f'(t) + r{t) - (2 + sin(2t) + cos(2t))f(t)^| 1^ ^ 0.1361 

Since it is enough for our purposes we can consider simpler rational approxima- 
tions of fo and fi, but keeping a similar accuracy. For finding these rational 
approximations, we search them doing the continuous fraction expansion of these 
values. For instance 

fo = [0,1,2,1,9,1,21,17,3,11] 

giving the convergents 1, 2/3, 3/4, 29/39, 32/43,. . .. Similarly fa gives 1/4, 1/5, 
28/139, 29/144,. . .. At this point we have the following new candidate to be an 
approximation of the periodic solution 

m = l-lcos{2t). 

Its accuracy w.r.t. equation (fT7|) is 



S= \ \f'{t)+f{t) - (2 + sin(2t) + cos(2t))f(t)^||, = ^ff^^^ ^ 0.1398 < 0.14, 



1600 



and so, quite similar to the one of r{t). 
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Therefore r{t) and f{t) are solutions of f fT7|) with similar accuracy so we keep 
f(t) as the second order approximation given by this modification of the HBM. 
For this f{t) we already know that its accuracy is 5* = 0.14. 

We need to know the value of M given in Theorem 11.11 With this aim we 
will apply Lemma (2.31 We consider in that lemma a function L{t) formed by 13 
straight lines and i = 1/9. Then we get that we can take M = 2.3. Therefore, 
since 2MS = 0.644 and 0.55 = § < f(t) < i = 0.95. 

We have that / = [-0.094, 1.594] in Theorem O Moreover 



i2 



< 6|2 + sin(2t) + cos(2t)||r| < (12 + 6V2)\r\ < y |r| 



Thus taking K = f(1.594) ^ 32.68 we get that 2M'^KS ^ 48.4 > 1 and we can 
not apply Theorem I l.li 

Doing similar computations with the successive approaches given by the HBM 
we obtain 

r{t) = I - |cos(2t) + ^cos(4t), 

f{t) = |-|cos(2t) + ^cos(4t)-^cos(6t). 
It is worth to comment that the above two functions are periodic functions that 
approximate to solution of (11 7p with accuracies 0.045 and 0.018, respectively, 
while the solutions obtained solving approximately the non-linear systems with ten 
significative digits have similar accuracies, namely 0.043 and 0.013, respectively. 
For none of both approaches Theorem 11.11 applies. Let us see that the next order 
HBM works for this example. 

If we do all the computations we obtain the candidate to be solution 

4 

f(t) = y^r2fcCos(2H), 

k=0 

with 



ro = 0.7457489122, ra = -0.2016836610, = 0.04065712547, 

re = - 0.009092599917, rg = 0.002133823488. 

Computing the accuracy of f(t) we obtain that it is 0.0039. If we take the approx- 
imation, using some convergents of r2k, 

m = 1-1 cos(2t) + ^ cos(4t) - ^ cos(6t) + ^ cos(8t) 

it has accuracy 0.0125. This means that we have lost significative digits and 
we need to take convergents of r2k that have at least 3 significative digits. For 
instance some convergents of ro are 1, 2/3, 3/4, 41/55, 44/59,. . . and we choose 
44/59. Finally we consider 

44 24 2 1 1 

'^'^ = 59 " 119 + 49 - no + 468 ^^^^ 

The accuracy of f is 0.00394 quite similar to the one of r{t). So we take S = 0.004. 
Let us see that Theorem 11.11 applies if we take this approximate periodic solution. 

In this case, by applying Lemma 12.31 using the piecewise linear function L 
formed by 10 pieces and £ = 1/10, we obtain that we can take M = 2.4. 
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Since it can be seen that 0.5 < r{t) < 1 and 2MS = 0.0192 we can take in 
Theorem O the interval / := [0.4808, 1.0192]. 
Then 



41 

< y(l-02) = 20.91 =: K. 



max 

/x[0,27r 

Finally, 2M'^KS ^ 0.96 < 1 and Theorem O applies 
Finally, it is easy to see that 

r.27r 



/ —X(f(t),t)dt> 12.5, 
Jo or 



which is bigger than 2ti/M 2.6. Therefore the hyperbolicity of the periodic 
orbit given by Theorem 11.11 follows. In short we have proved. 

Proposition 4.1. Consider the periodic function f(t) given in (12T|) . Then there 
is a periodic solution r*{t) of (|T71) . such that 

\\f-r*\\ < 0.0192, 

II I I oo — ' 

which is hyperbolic and unstable and it is the only periodic solution of ( I17p in this 
strip. 

Remark 4.2. Using the known analytic expression of r*(t) it can be seen that 
indeed 

\\f-r*\\ < 0.0007. 

II I I oo — 

Notice that by using a high enough HBM we have obtained a proof of the 
existence of a hyperbolic periodic orbit and an effective approximation f{t) without 
integrating the differential equation. 

4.2. A rigid cubic system. In this section we study some concrete cases of the 
family of rigid cubic systems 

X = -y - x{a + bx + x^), ^22) 
y = X — y{a + bx + x'^), ^ ' 

already considered in [2]. In that paper it is proved that (1221) has at most one 
limit cycle and when it exists is hyperbolic. With our point of view we will find 
an explicit approximation of the limit cycle, see Proposition 11.21 We consider the 
case a = —b = 1/10, that in polar coordinates writes as ([5]), 

dr 1 1 

/ /\2 2/\'"i 

^ ^ It ^ T6^ ~ To ^ ~ ^ ' 

and we start explaining how we have found the approximation of the periodic 
solution of ([5]) given in Proposition 11.21 

First attempt: the HBM. First we try to apply this method to find an approx- 
imation of the periodic solution of (|5]) that allows to use Theorem 11.11 

Searching for a solution of the form r(t) = tq and imposing that the first har- 
monic of 111 1 

-Tn rn H cos(t)rn H — cos(2t)rn 

2 ° 10 10 ^ ' ° 2 ^ ^ ° 

vanishes we obtain that 
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Hence tq = y/5/5 ~ 0.4472135954 is the first order solution given by the HBM. 
We obtain that the positive approximate solution is r = \/5/5. For applying the 
second order HBM we search for an approximation of the form 

r{t) = tq + ri cos(t) + Si sin(t). 

The vanishing of the coefficients of 1, cos(t) and sin(t) in J^{r{t)) provides the 
non-linear system 

9o 5o3 9 1q 3q lo 1 

-rnri ri H — risT H rr. H H s-, ri = 0, 

4 ° 8 ^ 8 ^ 10 ° 40 ^ 40 ^ 10 

-TnSi H — r. Si H — Si H ri Si Si — ri = 0, 

4 ° 8 ^ 8 ^ 20 10 

1q9o X 3o 1 

-rn H — rnri rn H — rnSi H rnri = 0. 

2 0^ g 1 10 ° 8 ° ^ 10 ° 

By using the same tools than in the previous example we obtain that one of 
the approximated solutions of the above system is ro ~ 0.4471066159, ri 
—0.0009814101 and Si ~ —0.0196567414. We search simple rational approxima- 
tions of ro, ri and Si, doing again the respective continuous fraction expansions 
and we obtain the candidate 

rW = ^-Y^cos(t)-lsin(t), 

to be an approximate periodic solution of (E]). It can be seen that it has accuracy 
S ~ 0.046. Doing all the computations needed to apply Theorem 11.11 we obtain 
that we are not under its hypotheses. Therefore we need to continue with the 
HBM of second order. 

Doing the second order approach we obtain five algebraic polynomial equations, 
that we omit for the sake of simplicity. Unfortunately, neither using the resultant 
method as in the previous cases, nor using the more sophisticated tool of Grobner 
basis, our computers are able to obtain an approximate solution to start our 
theoretical analysis. 

A numerical approach. First, we search a numerical solution of ([5]) by using the 
Taylor series method. From this approximation we compute, again numerically, 
its first Fourier terms, obtaining 

3 

f(t) = Tfc cos(H) + Sfc sin(fct), 

fc=0 

where 

ro = 0.4483561517, n = -0.0024133439, si = -0.0193837572, 
r2 = - 0.0037463296, S2 = -0.0220176517, 
rg = - 0.0012390886, sg = 0.0003784656. 

The accuracy of f(t) is 0.00289. If we take a new nicer approximation, using again 
some convergents of r^ and s^, we obtain 

r(t) = — cos(t) - — sin(t) — cos(2t) - — sin(2t) — cos(3t), (23) 

^ ^ 9 693 ^ ^ 51 ^ ^ 653 ^ ^ 45 ^ ^ 780 ^ ^' ^ ^ 

with accuracy 0.00298, quite similar to the one of fi(t). Note that fl23l) is precisely 

the approximation of the periodic solution of stated in Proposition 11.21 
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Proof of Proposition M.B. We already know that the accuracy of r(t) is 5* := 0.003. 
To apply Theorem II. II we will compute M and K. 
First we calculate A{t) = ^X{f{t),t). 

A (f\ _ 2891685439 _ 347888350813299559 f _ 561179 pf.o('+^ _ 685338551 • 
72733752000 1778094556332494400 36756720 ^''J 8000712720 

757058717 p^„|'94^ _ 40221206418131 ^\y,(')f\ _ 2923231 ^^^(Of-] 
48004276320 ^"-'=1^''; 273447836421760 ^^^H^f'J 576974475 

I 37724429 o^v.|'Q+^ _ 353400139 pf^^(Af\ i 17671001708653999 
36003207240 96008552640 '-^^l^'-J "r 42674269351979865600 ^^^'■K^'' J 

I 5358811 r^c^t^f] -L 4708003 ^i.^(^f] , 1537 ^^^(f,f] 
^300026727000 ^'^°\'-"'J ^ 20001781800 ^'■'''■y'-"' J ^ 207810720 ^'^^K^'^J 

+ li^^^'^oo sin(6t) + 327goo cos(7t) - -^^^^ sin(7t) 
~ 12979200 sin(8t). 

Now, by using again Lemma 12. 3^ we find a deformation constant M. In this 
case we use as lower bound for A the piecewise function L formed by 7 straight 
lines and i = 1/18. We obtain that we can take M = 7. Therefore 2MS ?a 0.042. 

Since it can be seen that 0.4 < f{t) < 0.47 in Theorem 11.11 we can consider the 
interval / = [0.358,0.512]. 

Then 

1 1 

<- + 6||f|L = - + 6(0.512) = 3.272 =: i^' 



max 

/x[0,27r] 



dr 



5 " "°° 5 



Finally, 2M'^KS ~ 0.962 < 1 and the first part of Theorem 11.11 applies. Hence 
equation ([5]) has a periodic solution r*(t) satisfying 

||r - r*|L < 0.042, (24) 



and is the only one in this strip. 
It can also be seen that 

^2n Q 







> 1.2. 



Since 2tt /M ^ 0.9, the hyperbolicity of r*(t) follows applying the second part of 
the theorem. □ 



Notice that the example of system f l22|) that we have studied is a = A and 6 = — A 
with A = 1/10. With the same techniques it can be seen that the same function 
f{t) given in the statement of Proposition 11.21 is an approximation of the unique 
periodic orbit of the system when |A — 1/10| < 1/500, which also satisfies 
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